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Alkali metals exhibit unexpected structures and electronic behavior at high pressures. Compres¬ 
sion of metallic sodium (Na) to 200 GPa leads to the stability of a wide-band-gap insulator with 
the double hexagonal hP4 structure. Post-hP4 structures remain unexplored, but they are impor¬ 
tant for addressing the question of the pressure at which Na reverts to a metal. Here we report 
the reentrant metallicity of Na at the very high pressure of 15.5 terapascal (TPa), predicted using 
first-principles structure searching simulations. Na is therefore insulating over the large pressure 
range of 0.2-15.5 TPa. Unusually, Na adopts an oP8 structure at pressures of 117-125 GPa, and 
the same oP8 structure at 1.75-15.5 TPa. Metallization of Na occurs on formation of a stable and 
striking body-centered cubic cI24 electride structure consisting of Nai 2 icosahedra, each housing at 
its center about one electron which is not associated with any Na ions. 

PACS numbers: 61.50.Ks, 61.50.Ah 62.50.-p, 71.20.-b 


Alkali metals have long been known to possess sim¬ 
ple electronic structures at ambient pressure that are 
well explained by a nearly-free-electron model. How¬ 
ever, the simple low-pressure structures of alkali met¬ 
als do not remain simple upon compression. Rich and 
complex phases and remarkable physical phenomena 
have been observed, such as greatly increased elec¬ 
trical resistivity [l]], enhanced superconductivity i], 
unusually low melting temperaturesM 5] and metal- 
insulator /semiconductor transitions [6J, 7 ]. 

The metal-insulator transitions in Na and lithium 
(Li) 00 are among the most fascinating observations 
at high pressure. Neaton and Ashcroft predicted that 
upon compression Li[s| and Na@] might transform into 
atomically-paired structures. If a paired structure was 
to be adopted, both Na and Li would have the potential 
to be semiconducting. These predictions have stimu¬ 
lated numerous studies of the high-pressure structures 
of Li and Na (e.g., Refs. @, 0, Er M |) ■ 

The metal-insulator transition in Na was established 
by a joint theoretical and experimental effort [j] and 
confirmed by further experiments [l 5| . Na is predicted 
and observed to transform into an optically transpar¬ 
ent phase at above 200 GPa[jj|. This phase was pre¬ 
dicted and experimentally confirmed to have a dou¬ 
ble hexagonal hP4 structure Q, with a remarkably 
large bandgap reaching 6.5 eV at 600 GPa. A metal- 
semiconductor transition was observed in Li at 80 GPa. 
The semiconducting structure of dense Li remained 
unsolved despite considerable attention (e.g., Refs. P 
0)- Blind crystal structure prediction calculations[12] 
on dense Li predicted a complex semiconducting oC40 
structure. Independent experimental and theoretical 
studies also found that semiconducting Li adopts the 
oC40 structure 0. 

Analysis of the charge density EH, 0 shows that 
the emergence of an insulating state in dense Na/Li is, 


however, not attributable to atomic pairing, but rather 
to strong localization of valence electrons within lattice 
voids. Notably, the insulating phases of Na and Li can 
be regarded as prototypical examples of electrides 17|, 
in which the ionic cores play the role of cations, and 
interstitial electrons form the anions 0, U, Li}. Stable 
electride structures have also been predicted in other 
elements at high pressures, e.g., A1 Q, Mg0, Ca0 
and K 2ll|. The formation of an electride in Na at 
high pressures is favorable because it reduces the ki¬ 
netic energy of the higher-energy electrons. These elec¬ 
tronic orbitals are forced to oscillate rapidly close to 
the atomic cores due to Pauli repulsion, but if they 
move away from the cores the oscillations and hence 
the kinetic energy can be reduced 17, HU . The elec¬ 
trides have open structures if one considers only the 
ions, but when the interstitial electrons are considered 
as anions the structures make good chemical sense as 
ionic solids. 

Pressure-induced metal-insulator transitions have 
been predicted in other materials at TPa pressures. 
For example, it was suggested that Mg, a good metal 
at normal pressures, would transform from a metal to 
a semimetal at 24 TPa[H|, and Ni was predicted to 
transform to an insulator at 34 TPaJHUj. The metal- 
insulator transition in Ni is driven by the complete 
4s—>-3d charge transfer under pressure, which is differ¬ 
ent from the electron localization found in Li and Na. 

Metallization is presumed to be the ultimate fate of 
all materials under sufficiently strong compression 241. 
The insulating phase of Ni has been predicted to re¬ 
vert to metallic behavior at 51 TPa[23|. Recent electri¬ 
cal resistance measurements have revealed that semi¬ 
conducting Li reverts to a ’’poor metal” above 120 
GPa0. The question of the pressure at which Na 
will also revert to a metal remains unanswered. De¬ 
termining the ground state structures of Na at high 
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pressures provides the key to whether insulating or 
metallic behavior is observed. At room temperature, 
Na undergoes pressure-induced phase transitions with 
an established sequence of bcc —> fee (65 GPa) [25j —> 
cI16 (103 GPa)[2<Il8| -s- oP 8 (117 GPa)@ tI19 
(125 GPa)[H -> hP4 (200 GPa) 0. The hP4 structure 
is the highest-pressure phase of Na known so far. 


In this letter, we report the reentrant metallicity of 
Na at a surprisingly high pressure of 15.5 TPa, pre¬ 
dicted via extensive structure searching in conjunc¬ 
tion with first-principles calculations. We predict that 
Na will show insulating behavior over a remarkably 
large pressure range of 0.2-15.5 TPa and will exhibit a 
phase transition to a second insulating orthorhombic 
oP 8 structure at 1.75 TPa. Metallization of Na ap¬ 
pears with the formation of a striking body-centered 
cubic (bcc) cI24 structure whose lattice sites are pop¬ 
ulated by Nai 2 icosahedra, each containing about one 
electron which is not associated with any Na ions. 


Our structure searching calculations with simulation 
cells containing up to 24 atoms were performed in a 
wide pressure range (0.4-20 TPah We used the efficient 
CALYPSO }29L i30| and AIRSS [3T], methods, which 

have both been successfully applied to investigating 
structures of materials at high pressures 12| [3l|, |33j- 
|42}. TPa pressures can now be achieved with dy¬ 
namic ramped compression. In particular, diamond 
has been studied experimentally up to 5 TPa[43|. The 
density functional theory[44j (DFT) calculations were 
performed using the CASTEP plane-wave code (45| 
with the Perdew-Burke-Ernzerhoffiij (PBE) general¬ 
ized gradient approximation (GGA) functional. We 
performed calculations with both norm-conserving and 
ultrasoft Na pseudopotentials, in which only the Is 
energy level was treated as core. An ultrasoft pseu¬ 
dopotential was used for the structure searches, with 
a plane-wave cutoff energy of 910 eV and a Brillouin- 
zone integration grid spacing of 27 r x 0.05A . The 

structures obtained were re-optimized at a higher level 
of accuracy for both the norm-conserving and ultra- 
soft pseudopotentials. Here, we present only norm- 
conserving results since the two pseudopotentials led 
to consistent conclusions. A cutoff energy of 1633 eV 
and a k-point grid spacing of 2n x 0.03A were used 
in the calculations with norm-conserving pseudopoten¬ 
tials. The phonon dispersion curves were computed by 
a finite displacement method as implemented in the 
CASTEP code. 


Our structure searches readily reproduced the 
known double hexagonal hP4 structure (space group 
P 63 /mmc, 4 atoms/unit cell) above 400 GPa and we 
found no better structure below 1.5 TPa. Surpris¬ 
ingly, above 2 TPa, we found the reappearance of 
the orthorhombic 0 P 8 structure (space group Pnma, 
8 atoms/unit cell, Fig. [lji) which is also observed at 
low pressures of 117-125 GPa[26|. Strikingly, we pre¬ 
dict a hitherto unexpected highly symmetric bcc struc¬ 
ture at 20 TPa (space group Im-3, 24 atoms/unit cell, 
denoted cI24 hereafter, see Fig. lb). The cI24 struc¬ 
ture consists of Nai 2 cages, each of which contains 12 



FIG. 1: (color online) Crystal structures of 0 P 8 (a) and 
cI24 (b) and their electron density difference maps (c-d) at 
ultra-high pressures. The lattice parameters of oP 8 -Na at 
2 TPa are a = 2.88 A, b = 2.12 A, and c = 4.0 A with 
two inequivalent Nal and Na2 atoms at 4c positions of 
(0.809, 0.75, 0.567) and (0.523, 0.25, 0.719), respectively. 
For c.I24-Na at 16 TPa, a = 3.12 A with Na atoms at 24g 
positions (0, 0.302, 0.183). Electron density differences are 
plotted within the (010) plane of oP 8 -Na (c) at 2 TPa and 
cI24-Na (d) at 16 TPa, respectively. Al and A2 in (d) 
represent electron attractors located at the centers of the 
Nai 2 icosahedra and octahedral voids formed by the six 
neighboring Nai 2 icosahedra, respectively. 


Na atoms forming an icosahedron. The shortest intra- 
and inter-icosahedral Na-Na distances at 16 TPa are 
calculated to be 1.138 A and 1.168 A, respectively, in¬ 
dicating fairly strong core-core overlap (the ionic ra¬ 
dius of Na + is 1.02 A). An analysis of electron density 
differences (Fig. QJl) reveals that the cI24 phase is also 
an electride with two distinct electron attractors lo¬ 
cated at the centers of the Nai 2 icosahedra (Al) and 
the octahedral voids (A2). However, unlike the nearly 
spherical isosurface of the interstitial electrons in the 
hP4 and 0 P 8 structures, dumbbell-like electron local¬ 
ization is found in the octahedral voids. 

The enthalpies of the 0 P 8 (relative to hP4) and cI24 
(relative to 0 P 8 ) phases are plotted as functions of 
pressure in Figs. [2] (a) and (b), respectively. The in¬ 
set to Fig. [2] (a) shows that hP4-Na transforms to 
the 0 P 8 structure at 1.75 TPa with a continuous vol¬ 
ume change, indicating a second-order phase transi¬ 
tion. The 0 P 8 phase remains stable up to 15.5 TPa, 
above which the cI24 phase becomes the most stable. 
Phonon calculations for 0 P 8 and cI24-Na in their cor¬ 
responding stable pressures ranges show no imaginary 
frequencies, demonstrating the dynamical stability of 
the structures (Fig. [5£, d). The maximum frequency 












3 




W 0-2 



hP4 




\ 1 


oP8 


\ ' N. 

Q 0.0 

X 

^ 4 

\ ops 


N. 

E 3 

.hp+Xj'™ \ 



£ 2 

: 

*- 0.2 

<1 

cI24 


1 Pressure (T$a) ^ ■ 


♦ 




[OllOJ 


( d ) 



FIG. 2: (color online) Enthalpy differences (relative to hP4 
(a) and oP8 (b)) of Na as a function of pressure. The inset 
in (a) shows the pressure dependence of the volumes of the 
hP4 and oP8 phases. The predicted sequence of phase tran¬ 
sitions at zero temperature in the pressure range of 0 - 20 
TPa is provided in supplementary Figure SI. The volume 
appears to change continuously across the phase transition, 
(c) and (d) show the calculated phonon dispersion relations 
of oP8-Na at 2 TPa and cI24-Na at 16 TPa, respectively. 


in the TPa pressure regime reaches 137 THz at 16 
TPa, which is very high and implies a large nuclear 
zero-point (ZP) energy. The estimated ZP energies 
within the quasi-harmonic approximation are indeed 
extremely large: 0.577 and 0.585 eV/atom for cI24 and 
oP8 phases at 20 TPa, respectively. However, the dif¬ 
ference in ZP energies of only 8 meV/atom is too small 
to modify the phase diagram of Na at TPa pressures. 
Furthermore, we have calculated the Gibbs free ener¬ 
gies of the cI24 and oP8 phases at 20 TPa within the 
quasi-harmonic approximation at finite temperatures 
to account for vibrational contributions. The resul¬ 
tant difference in free energy between cI24 and oP8 at 
1000 K is only 3 meV/atom larger than that at 0 K 
(Fig. (2b), which gives a negligible change in the oP8 - 
cI24 transition pressure. 

We calculated the phonon dispersion curves of hP4- 
Na at various pressures in order to understand the 
origin of the re-emergence of the oP8 phase and the 
second-order nature of the hP4 - oP8 transition. We 
found that a transverse acoustic (TA) phonon mode 
at the zone boundary M (0 0.5 0) point softens with 
increasing pressure (Fig. 3a) and its frequency goes to 
zero at 1.72 TPa (Fig. 3b). A frozen-phonon calcu¬ 
lation was performed by distorting the hP4 structure 
according to the atomic vibrations for the TA mode at 
the M point as depicted in Fig. EIc), and structural 
optimizations were subsequently carried out. The re¬ 
sultant energetically stable structure is indeed the or¬ 
thorhombic oP8 phase, as indicated by the dashed cell 
in Fig. [3] (d). It is clear that the hP4 - oP8 tran¬ 
sition is driven by the soft TA mode. This result is 
consistent with the second-order hP4 - oP8 transition 
characterized by a continuous volume change (inset in 
Fig. (2^,). Soft-phonon driven hP4 - oP8 transforma- 


FIG. 3: (color online) Phonon dispersion curves of hP4-Na 
at 2 TPa (a) and the phonon frequencies of the softened 
TA phonon mode at M (0, 0.5, 0) as a function of pressure 
(b). (c) shows the atomic vibrations of the soft TA mode at 
the M point viewed along the a-axis in a 1x2x1 supercell 
of the hP4 structure. Here, the Nal (Na2) and Na4 (Na3) 
atoms vibrate parallel (antiparallcl) to the [01 10] direction, 
while the Na5 (Na6) atoms antiparallel (parallel) to the 
[0001] direction. The gray and black balls indicate atoms 
in different planes. Arrows indicate vibrational direction 
of different Na atoms as labeled and are plotted from the 
eigenvector of the TA soft mode at the M point. The blue 
dashed cell in (d) indicates the oP8 structure resulting from 
the distortion. 



FIG. 4: (color online) Electronic band structure and pro¬ 
jected density of states of oP8-Na at 1.75 (PDOS, in units 
of eV -1 per atom) (a) and cI24-Na at 16 TPa (b). The 
horizontal dotted lines indicate the Fermi level, (c) and 
(d) show the variation of band gaps with pressure for hP4 
and oP8-Na, respectively, using the HSE06 functional. 


tions have also been reported in MnAs, CrTe, MnTe, 
MgTe, CrSb and CaH 2 • 

The electronic band structures and band gaps of 
hP4, oP8, and cI24-Na were calculated at selected 
pressures using the HSE06 functional IHH, which 
is known to give better band gaps than the standard 
PBE functional. The band gap of hP4-Na increases 
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with pressure and reaches 9.5 eV at the transition to 
the 0 P 8 phase at 1.75 TPa. At the transition, oP 8 -Na 
is also a wide band gap insulator with an even larger 
gap of 9.7 eV (Fig. [Hr). The gap of oP 8 -Na decreases 
with pressure down to 1.6 eV at 15.5 TPa. The oP 8 
structure is compressed by a factor of three from 125 
GPa (volume = 77.9 A 3 per unit cell) to 1.75 TPa (vol¬ 
ume = 25.9 A 3 per unit cell). As a result of this strong 
compression, the oP 8 structure at TPa pressures forms 
an electride with valence electrons trapped at the in¬ 
terstitial sites (Fig. QJ:). The electronic structure is 
fundamentally different from that of the nearly-free- 
electron behavior of oP 8 at a much lower pressure of 
125 GPa (supplementary Figure S2)[j52|]. As a result, 
the oP 8 phase is insulating at TPa pressures, but is 
metallic at pressure below 125 GPaI26l I52|. 

After the transition to the cI24 phase, the localized 
electrons in the octahedral voids are partially squeezed 
out due to the decreasing interstitial space, forming a 
series of conducting channels between A 2 attractors 
and Na cations, as shown in Fig. |T| (d). Notably, 
the band structure of cI24-Na at 16 TPa shows highly 
dispersive bands crossing the Fermi level (Fig. [3 )d), 
and the return to metallicity. 

The enthalpy H = U+PV governs the phase sta¬ 
bility at 0 K. Under sufficiently strong compression 
the hard core repulsion between atoms becomes very 
important and it favors close-packed structures such 
as (bcc, fee or hep) @ which have small volumes V. 
However, the cI24 structure which Na adopts at TPa 
pressure is far from simple, and is constructed from 
Nai 2 icosahedral cages, further challenging the estab¬ 
lished picture of simple packing at high densities. In¬ 
deed, our optimized hypothetical bcc structure with 
Na atoms occupying the lattice sites at 20 TPa gives 
a volume of 1.175 A 3 /atom, which has a lower den¬ 
sity than cI24-Na (1.148 A 3 /atom) whose bcc lattice 
sites are occupied by Nai 2 icosahedral cages. Integra¬ 
tion of the electron density at 16 TPa shows that the 
charge on each Na atom is depleted by about ^0.35 
electrons, while the cage-center Al and interstitial A2 
sites (Fig. [TJl) have excess electronic charges of about 
~1.2 and ~1.0 electrons, respectively. As a whole, 9 of 
the 12 valence electrons remain within each Nai 2 cage 
(including the localized central electrons), while 3 elec¬ 
trons localize in the octahedral voids. The 9-electron 
configuration might allow an Nai 2 cage to be termed a 
super-alkali atom. Electronically, a Nai 2 cage resem¬ 
bles a Potassium atom with an electronic configuration 
of 3s 2 3 p e 4s 1 . 

Cage-like structures are seldom seen at high pres¬ 
sures since they do not usually favor dense packing. 
Known cage structures are typically found in elements 
in which it is easy to form covalent bonds. Exam¬ 
ples occur in various allotropes of B (a-B^jHiJ and 7 - 
B 2 8 551) and C (Ceo and C 70 ) [56]. and Si/Ge cages in 


clathrates (e.g., in Ba8Si46[57|). In recent work we also 
found an N 10 cage structure in polymeric nitrogen[58j] 
and an H 2 4 cage in hydrogen-rich CaH fi [37]. Here, 
the formation of Nai 2 icosahedra is rather surprising 
since Na is unable to stabilize covalent bonds. The 


Nai 2 cage is geometrically similar to a B 12 cage, which 
is the fundamental building block of various boron 
allotropes [mI [55| . B is electron deficient and forms 
multi-center covalent bonds. Three-center covalent 
bonds within the icosahedra and two/three-center co¬ 
valent bonds between the icosahedra in various struc¬ 
tures of elemental B satisfy the octet rule and generate 
insulating states. Here, there is no covalent bonding 
within the Nai 2 cages and it is actually a well-packed 
ionic structure. The electrons at the center of the 
cage contribute to the stability of the Nai 2 cages, and 
play a similar role to that of Ba atoms in the BasSUe 
clathrate structure [571. 

We find that cI24-Na is structurally similar to the 
bcc form of the Ali 2 W alloy [mJ, in which AI 12 also 
forms an icosahedral cage with a W atom localized at 
its center, and the W-encaged Al 12 icosahedrons oc¬ 
cupy bcc lattice sites. One AI 12 W unit resembles a 
Nai 2 icosahedron, and the encaged electron plays the 
role of the pseudo-ions in W. In contrast to the single 
dumbbell-like electron attractor found in the octahe¬ 
dral voids of cI24-Na, the voids in AI 12 W contain four 
distinct dumbbell-like electron attractors (supplemen¬ 
tary Fig. S5)[H. 


Na is insulating over the large pressure range of 0.2- 
15.5 TPa. The pressure interval of 15.3 TPa (98.7% of 
the pressure range 0 - 15.5 TPa) is comparable to the 
17 TPa (33.3% of the whole pressure range) pressure 
interval predicted for the insulating phase of Nij23i]. 
The pressure for re-entrant metallicity in Na (15.5 
TPa) is much higher than in Li (120 GPa), which may 
be associated with the entirely different mechanism for 
insulating behavior. In the hP4 and 0 P 8 structures 
of Na, the valence electrons localize strongly within 
the lattice interstices and maintain a nearly spherical 
shape (Fig. [lj:) up to 15.5 TPa. However, semicon¬ 
ducting Li has a complex oC40 structure with intersti¬ 
tial regions of different shapes in which electrons are 
localized at a rather low pressure of 80 GPa. These lo¬ 
calized electrons are connected to each other and form 
conducting channels at 120 GPa (52]. 
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